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Abstract 

Cramer's rules for some left, right and two-sided quaternion ma- 
trix equations are obtained within the framework of the theory of the 
column and row determinants. 
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1 Introduction 

Quaternion matrix equations play important roles in both theoretical studies 
and numerical computations of quaternion application disciplines [HEIIH] and 
have been studied by many experts [5l [6l [?j in mainly by means of complex 
representation of quaternion matrices. In this paper we obtain Cramer's 
rule for some left, right and two-sided quaternion matrix equations within 
the framework of the theory of new matrix functionals over the quaternion 
skew field (the column and row determinants) introduced in [1]. In the first 
point we shortly cite some provisions from the theory of the column and row 
determinants which are necessary for the following. The theory of the column 
and row determinants of a quaternion matrix are considered completely in [1]. 
In the second point the Cramer rules for left, right and two-sided quaternion 
matrix equations are obtained. 
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2 Elements of the theory of the column and 
row determinants. 

Let M (n, H) be the ring of n x tt, quaternion matrices. By H™^" denote 
the set of all m x n matrices over the quaternion skew field H and by H^^" 
denote its subset of matrices of rank r. Suppose Sn is the symmetric group 
on the set /„ = {!,..., n}. 

Definition 2.1 The ith row determinant of A = {uij) G M (n, H) is defined 
by 



rdetjA — ^ ( 1) ajj^^ai^,^i^,^^^ . . .aj^^^^^i . . . a, 
o-eSn 



^kr^kr+l ■ ■ ■ ^hr+lr^kr 



for all i = 1, . . . ,n. The elements of the permutation a are indices of each 
monomial. The left- ordered cycle notation of the permutation a is written as 
follows, 

The index i opens the first cycle from the left and other cycles satisfy the 
following conditions, ik^ < ik^ < ■ ■ ■ < hr ^'<^d i^^ < ikt+s for all t — 2, . . . ,r 
and s — 1, . . . ,lt. 

Definition 2.2 The jth column determinant of A = (aij) e M{n,M) is 
defined by 

cdetj A = (~1) %fc^j)c^+i^ • • • %fc^+ii)c^ • • ■O'jjk^+i^ ■ ■ ■ O'jk^+ijk^CLjk^j 

for all j — l,...,n. The right-ordered cycle notation of the permutation 
T & Sn is written as follows, 

'T = iJkr+lr ■ ■ -jkr + dkr) ■ ■ ■ {jki+h ■ ■ ■ j ki+lj ki) Uki+h ■ ■ -jki+ljkij) ■ 

The index j opens the first cycle from the right and other cycles satisfy the 
following conditions, jk^ < jka < ■ ■ ■ < jkr o^^d j^t < jkt+s for all t = 2, . . . ,r 
and s — 1, . . . ,lf 
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Suppose A*-' denotes the submatrix of A obtained by deleting both the 
ith row and the jth column. Let a j be the jth column and aj be the ith 
row of A. Suppose A (b) denotes the matrix obtained from A by replacing 
its jth column with the column b, and Aj. (b) denotes the matrix obtained 
from A by replacing its ith row with the row b. 

We note some properties of column and row determinants of a quaternion 
matrix A = (ajj), where i e In, j G J„ and /„ = J„ = {1, . . . , n}. 

Proposition 2.1 If b E M, then rdetjAj, (6- a^.) = b ■ rdetjA for all 
i = 1, . . . ,n. 

Proposition 2.2 If b E M, then cdetj A j (a jb) = cdetj A6 for all j = 
1, . . . ,n. 

Proposition 2.3 If for A G M (n, H) there exists t E In such that atj = 
bj + Cj for all j = 1, . . . ,n, then 

rdetj A = rdetj At, (b) + rdet. At, (c) , 
cdetj A = cdetj A^, (b) + cdetj A^. (c) , 

where b = . . . , c = (ci, . . . , c„) and for all i = 1, . . . ,n. 

Proposition 2.4 If for A E M (n, H) there exists t E Jn such that 
(lit = bi + Ci for all i = 1, . . . ,n, then 

rdetj A = rdet j A t (b) + rdetj A ^ (c) , 
cdetj A = cdetj A t (b) + cdetj A . t (c) , 

where b = (6i, . . . , bn)^ , c = (ci, . . . , Cn)^ and for all j = 1, . . . , n. 

The following lemmas enable us to expand rdetj A by cofactors along the 
ith row and cdetj A along the jth column respectively for alH, j = 1, . . . , n. 

Lemma 2.1 |^ Let Rij be the right ij-th cofactor of A eM {n,M), that is, 

n 

rdetj A = ^ Ojj ■ Rij for all i = 1, . . . ,n. Then 

^ f -rdetj A^i (a.,) , z ^ j, 
\ rdetfcA", i = j, 

where A" (a. j) is obtained from A by replacing the jth column with the ith 
column, and then by deleting both the ith row and column, k = min {/„ \ {i}}. 
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Lemma 2.2 JJJ Let Lij be the left ij-th cofactor o/ A G M{n,M), that is, 

n 

cdetj A = ^ Lij ■ ttij for all j = 1, . . . ,n. Then 
1=1 

T ^ / -cdeti A[^ (aj.) , j, 
\ cdetfcA"-, 1 = 3, 

•where A]-' (aj.) is obtained from A by replacing the ith row with the jth row, 
and then by deleting both the jth row and column, k = mm{J„ \ {j}}. 

We recall some well-known definitions. The conjugate of a quaternion a — 
Uq + aii + a2j + 03 A; G EI is defined by a = Oo — oiz — 02 j — a^k. The Hermitian 
adjoint matrix of A = (ajj) G H"^™' is called the matrix A* = if 
a*j = CLji for alH = 1, . . . , n and j = 1, . . . ,m. The matrix A = (aij) G H"^*" 
is Hermitian if A* = A. 

A following theorem has a key value in the theory of the column and row 
determinants. 

Theorem 2.1 If A = (uij) G M(n,H[) is Hermitian, then rdetiA = 

■ ■ ■ = rdet„A = cdetiA = ■ ■ ■ = cdet„A G M. 

Taking into account Theorem 12.11 we define the determinant of a Hermitian 
matrix by putting det A := rdetjA = cdet^A for all i = l,...,n. This 
determinant of a Hermitian matrix coincides with the Moore determinant. 
The properties of the determinant of a Hermitian matrix are considered in 
[5] by means of the column and row determinants. Among them we note the 
followings. 

Theorem 2.2 ^ // the ith row of a Hermitian matrix A G M{n,M.) is 
replaced with a left linear combination of its other rows, i.e. a^. = ciaj^. + 

■ ■ ■ + Cfcajj.., where q G H for all I = 1, . . . , k and {i, ii} C /„, then 

cdetiAj. (ci ■ ajj. H h ■ a^^. ) = rdetjAj. (ci ■ a^,. H h ■ a^^. ) = 0. 

Theorem 2.3 //^/ // the jth column of a Hermitian matrix A G M {n,M.) 
is replaced with a right linear combination of its other columns, i.e. a^ = 
a jjCi + ■ ■ ■ + Si ji^Ck, where q G EI for all I = 1, . . . , k and {j,ji} C J„, then 

cdetiA. j (a.ij ■ ci H h a.^^ ■ c^) = rdetjA.j (a.j^ ■ ci H h a.^^ ■ c^) = 0. 
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The following theorem about determinantal representation of an inverse ma- 
trix of Hermitian follows immediately from these properties. 



Theorem 2.4 |^/ // a Hermitian matrix A G M (n, H) is such that det A ^ 
0, then there exist a unique right inverse matrix {RA)^^ and a unique left 
inverse matrix {LA)~^, and {RA)^^ = (LA)^'^ =: A^^. They possess the 
following determinantal representations : 



(RA)-' 



det A 



(LA) 



-1 



1 



det A 



/Ru 


R2I ■ 


Rnl^ 


Rl2 


R22 ■ 


Rn2 


\Rln 


R2n ■ 


Rnn ) 


fLn 


L2I ■ ■ 






L22 ■ ■ 


Ln2 


\Lln 


L2n ■ ■ 


Lnn ) 



(1) 



(2) 



Here Rij, Lij are right and left ij-th cof actors of A respectively for all i,j = 
1, . . . ,n. 

To obtain determinantal representation of an arbitrary inverse matrix A~^, 
we consider the right AA* and left A*A corresponding Hermitian matrix. 

Theorem 2.5 If an arbitrary column of A E M"*^" is a right linear 
combination of its other columns, or an arbitrary row of A* is a left linear 
combination of its others, then det A*A = 0. 



Since the principal submatrices of a Hermitian matrix are Hermitian, the 
principal minor may be defined as the determinant of its principal submatrix 
by analogy to the commutative case. We introduce the rank by principal 
minors that is the maximal order of a nonzero principal minor of a Hermitian 
matrix. The following theorem determines a relationship between it and the 
rank of a matrix defining as ceiling amount of right- linearly independent 
columns or left-linearly independent rows which form basis. 

Theorem 2.6 //^/ A rank by principal minors of A* A is equal to its rank 
and a rank of A e H™-^". 
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Theorem 2.7 |^ // A G H"^^", then an arbitrary column of A is a right 
linear combination of its basis columns or an arbitrary row of A is a left 
linear combination of its basis rows. 

The criterion of singularity of a Hermitian matrix is obtained. 

Theorem 2.8 //^/ The right- linearly independence of columns of A & jgjmxn 
or the left-linearly independence of rows of A* is the necessary and sufficient 
condition for det A* A ^ 0. 

Theorem 2.9 /^Z // A e M {n, H), then det AA* = det A*A. 

A concept of the double determinant is introduced by this theorem. This 
concept was initially introduced by L. Chen in [2]. 

Definition 2.3 The determinant of the corresponding Hermitian matrix of 
A G M(n,EI) is called its double determinant, i.e. ddetA := det (A* A) = 
det(AA*). 

The relationship between the double determinant and the noncommutative 
determinants of E. Moore, E. Study and J. Diedonne is obtained, ddetA = 
Mdet (A* A) = Sdet A = Ddet^A. But unlike those, the double determinant 
can be expanded along an arbitrary row or column by means of the column 
and row determinants. 

Definition 2.4 Suppose A G M(n,EI). We have a column expansion of 
ddetA along the jth column, ddetA = cdetj (A* A) = ' (^ij> '^'^^ ^ ''^^w 



expansion of it along the ith row, ddetA = rdetj (AA*) = ^ ajj-Mjj for all 

j 

i, j = 1, . . . ,n. Then by definition of the left double ijth cofactor we put hij 
and by definition of the right double ijth cofactor we put Mj^. 

Theorem 2.10 [4^ The necessary and sufficient condition of invertibility of 



A G M(n,e) IS ddetA ^ 0. Then there exists A'^ = (LA)~^ = (i?A)"\ 



where 



{LA) 



1 



(A* A)"' A* 



1 




(3) 



ddetA 
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(«A)-' = A-(AAT' = 5dSA 



and hij = cdetj(A*A).j (a*j), Mij = rdetj(AA*)i. (a*) for all i,j 





^ 11 


R 21 . . 




i 


^12 


K 22 • • 


K ri2 




^ In 


R 2n • • 


nnj 



(4) 



n. 



This theorem introduces the determinantal representations of an inverse ma- 
trix by the left ([1]) and right (|2]) double cofactors. The inverse matrix A~^ of 
A G M(n,EI) on the assumption of ddetA 7^ is represented by the analog 
of the classical adjoint matrix. If we denote this analog of the adjoint matrix 
over H by Adj[[A]], where Adj[[A]] = (Lij)nxn or Adj[[A]] = (Mij)nxn, then 
the following formula is valid over H: 



A 



-1 Adj[[A]] 



ddetA 



Using the determinantal representations of an inverse matrix by the left ([T]) 
and right ([2]) analogs of a classical adjoint matrix we obtain the Cramer rule 
for right and left systems of linear equations respectively. 

Theorem 2.11 Let A-x = y be a right system of linear equations with a ma- 
trix of coefficients A G M(n, H), a column of constants y = (yi, . . . , y„)^ e 
jjnxi ^ column of unknowns x = (xi, . . . , Xn)^ ■ If ddetA 7^ 0, then we 
have for allj = l,n 

cdetj(A*A).j(f) 



ddetA 



where f = A*y. 



Theorem 2.12 Let x • A = y &e a left system of linear equations with a 
matrix of coefficients A G M{n,M.), a row of constants y = {yi, . . . ,yn) G 
jgjixra ^ of unknowns x = (xi, . . . ,x„). //ddetA 7^ 0, then we obtain 
for all i = l,n 



rdet,(AA*)i. (z) 



ddetA 
where z = yA*. 
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3 Cramer's rule for some matrix equations. 



We denote A*B =: B = (bij), BA* =: B = (kj). 
Theorem 3.1 Suppose 

AX = B (5) 

is a right matrix equation, where {A,B} G M{n,Il) are given, X G M(n,M.) 
is unknown. If ddet A ^ 0, then ^ has a unique solution, and the solution 
is 

cdeti(A*A).i fb.,) 

x^, = (6) 

ddetA 

where h j is the jth column o/B for all i,j = 1, ...,n. 

Proof. By Theorem 12. 101 the matrix A is invertible. There exists the unique 
inverse matrix A~^. From this it foUows that the solution of (|5]) exists and 
is unique, X = A^^B. If we represent A~^ = (A*A)~^A* as a left inverse 
by ([3]), and use the determinantal representation of (A*A)~^ by ([2]), then for 
all i,j = 1, n we obtain 

1 " 

k=l 

where Lij is a left ijth cofactor of (A* A) for all i, j = 1, t?,. From this by 
Lemma [2.21 and denoting the j-th column of B by b it follows (^.M 

Theorem 3.2 Suppose 

XA = B (7) 

is a left matrix equation, where {A,B} G M(n,EI) are given, X G M{n,M.) 
is unknown. //ddetA 7^ 0, then ^ has a unique solution, and the solution 
is 

rdet,(AA*),-. (b,.) 
ddetA 

where hi is the ith column o/B for all i,j = 1, ...,n. 

Proof. By Theorem 12.101 the matrix A is invertible. There exists the unique 
inverse matrix A~^. From this it follows that the solution of ([7]) exists and is 
unique, X = BA~^. If we represent (A)^^ = A* (AA*)~^ as a right inverse 
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by (jl]) and use the determinantal representation of (AA*) ^ by ([T]), then for 
all i,j = 1, n we have 

1 " 

k=l 

where Rij is a right ijth cofactor of (AA*) for all i,j = 1, ...,n. From this 
by means of Lemma [2.11 and denoting the zth row of B by bj ., it follows ()8]). 
■ 

We denote A*CB* =: C = (5^). 
Theorem 3.3 Suppose 

AXE = C (9) 

is a two-sided matrix equation, where {A,B,C} G M(?t,, H) are given, X G 
M(n,E[) is unknown. //ddetA ^ and ddetB 7^ , then ^ has a unique 
solution, and the solution is 

_ rdetj(BB*),-. (c^) 
~ ddetA -ddetB ' 

or 

_ cdet,(A*A)., (eg) 

ddetA -ddetB ' ^ ^ 

where cf- := (cdetj(A*A). j (c.i) , . . . , cdetj(A*A). j (c.„)) is the row vector 
and := (rdetj(BB*)j. (c^^ ) , . . . , rdetj(BB*)j. (c„ ))"^ is the column vector 

and Ci,, c.j are the ith row vector and the jth column vector ofC, respectively, 
for all i,j = 1, n. 

Proof. By Theorem 12 . 1 01 the matrices A and B are invertible. There exist the 
unique inverse matrices A~^ and B~^. From this it follows that the solution of 
(jg]) exists and is unique, X = A"^CB"\ If we represent A"^ = (A*A)"^A* 
as a left inverse and (B)^^ = B*(BB*)~^ as a right inverse, then for all 
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i,j = 1, ...,n we have 



X 



X 



/ Xii Xi2 
^21 X22 



/ Cll C12 

C2I C22 



\Cnl C„2 



(A*A)-iA*CB*(BB*'"^ 



/ jA J A 



1 

ddctA 



Cln \ 
C-nn j 



tA tA 

-^12 -^22 



1 



A T A 

In '^2n 

-^12 ^22 



ddetA 



\R^n R 



2n 



^n2 

LnnJ 
R^l\ 

^n2 

dB / 



where is a left zjth cofactor of (A* A) and i?^- is a right zjth cofactor of 
(BB*) for alH, J = 1, n. This imphes 



m=l \k=l 



R 



ddetA ■ ddetB 



(12) 



for all i, j = l,n. From this by Lemma \2.2\ we obtain 

n 

^Lf-Ckm = cdeti(A*A).i (c.^) , 



k=l 



where c is the mth column- vector of C for all m = 1, ...,n. Denote by 
:= (cdeti(A*A).j (c.i) , . . . , cdeti(A*A). j (£.„)) the row-vector for all i = 



1, ...,n. Reducing the sum ^ ( ^ L^-Ckm ]R^m Lemma [TTl we obtain 

m=l Va;=1 / 

an analog of Cramer's rule for ([9]) by ffTOj) . 

Having changed the order of summation in f|T2l) . we have 



k=l 



m=l 



'jm 



ddetA ■ ddetB 



By Lemma [2?T| we obtain ^ CkmRfm = rdetj(BB*)j. (c^ ,), where c^, is 



m=l 

a kth row- vector of C for all A; = 1, n. We denote by 



(rdetj(BB*)j. (ci.) , . . . , rdetj(BB*)j. 
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the column- vector for all J = l,...,n. Reducing the sum ^ L^- ( ^ CkmRfm 

k=l \m=l J 

by Lemma [2.2[ we obtain Cramer's rule for ([9]) by f|TT]) . ■ 

In solving the matrix equations by Cramer's rules (jS]), dH]), fflU]) - f lTT]) we 
do not use the complex representation of quaternion matrices and work only 
in the quaternion skew field. 



4 Example 

Let us consider the two-sided matrix equation 

AXB = C (13) 

-j k 







j 


A 






i 






i 


k 


i 


1 and C = 1 


k 


j 








1 








1 
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6 j + 3k -j - fc^ 
we have A* = \ j i -k \ , A* A = \ -j - 3k 3 i | and 

j + k — z 3 

-i j \ ( ^ — ^ + 

-i -k 1 , BB* = 3i - J 3 l + 2k 
J I -2k 3 

2k 1-i-k 3 + 2 + 3A;\ 
-2 — Ai —2 + i ~ k i — k .It is easy to get, 
-4 + 2i l + 2i + j 1 + 4i + j J 
ddetA = det A* A = 8 and ddetB = det BB* = 4. Therefore has a 
solution. We shall find it by (JTOl) . At first we obtain the row- vectors cf- for 
all i = 1,2,3. 



/ 2k j + 3k -j-k\ 
cdeti(A*A).i(c.i)=cdeti -2-4i 3 i = 3 ■ 3(2A;)- 

\-A + 2i -i 3 ) 

-i{-i){3j + 5/e) + (-i - k){-i){-2 - M) - 3{j + 3A;)(-2 - 4z) + 
+ (i + 3k)i{-4: + 2i) - 3{-j - A;) (-4 + 2t) = 
= 24j + 8k, 
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and so forth. Continuing in the same way, we get 

cf, = (24j + 8A;, -8 - 8i + Aj + Ak,8 + 8i + Aj + Ak), 
= (-20 - 36i, -10 ~2i- 12j - 12k, -2 - 2i + 12j + Ak), 
cf = (12 + Ai,Q + 2i + 12j -Ak,Q + lOi - Aj + Ak). 

Then by ( [T0|) we have 

rdcti(BB*)i. (cf ) 
■^11 ~ ddctA-ddetB ~ 

/2Aj + 8k -8-8i + Aj + Ak 8 + 8i + Aj + Ak 
= ^ ■ rdeti I 3i-j 3 1 + 2A; 

\ i-j I -2k 3 

= ^ ■ {{2Aj + 8A;) • 3 ■ 3 - (24j + 8A;)(1 + 2A;)(1 - 2k) + 
+ (-8 -8i + Aj + 4A;)(1 + 2k){i - j) - (-8 - 8i + Aj + 4A;)(3z - j)3+ 
+ (8 + 8i + Aj + 4A;)(1 - 2A:)(32 - j) - (8 + 8z + Aj + Ak){i - i)3) = 

= ^ ■ (-32 + 32i), 

and so forth. Continuing in the same way, we obtain 

^ -32+32i ^ -88-72i+24i-8fc ^ 24+8i-40j+56fc 

11 32 ' *-^12 22 ' "^13 32 ' 

_ -16i+32j-48fc _ 20-28t-116j-76fc _ -44+68i+20j+12fc 

21 32 7 ^^22 32 5 -*'23 32 ? 

_ 16+16j+32fc _ 20+44i+52j-28fc _ -12-20i+12j-4fc 

'^l ~ 32 ' 32 — 32 ) -t'32 — 32 
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